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Abstract 



The nonlinear dynamics of dissipative quantum systems in incoherent laser fields is 
studied in the framework of master equation with random telegraph model describing the 
laser noise and Markovian approximation for dealing with the system-bath couplings. Floquet 
theory and time-dependent perturbation methods are used to facilitate both analytical and 
numerical solutions. We develop a theoretical formalism that provides a powerful tool for the 
detailed analysis of the dissipative quantum dynamics of multilevel systems driven by intense 
stochastic laser fields. It is found that the system relaxes to a steady state by the effect of laser 
phase and frequency noises and the kinetic of this relaxation increases with the addition of 
dissipative terms, introduced by the coupling to the reservoir. The case of amplitude 
fluctuations shows a different behavior. Other results concerning the destruction of quantum 
coherence and the dynamical localization will be established and further relaxation 
mechanisms such spontaneous emission and ionization process will be also considered. 



I- INTRODUCTION 



The study of the time evolution of quantum systems, which are on the one hand driven 
by an external field and in contact with a heat bath (reservoir) on the other hand, has received 
a great deal of attention in recent years [1-4]. In quantum optics, such systems are 
investigated in the dressed-atom picture of resonance fluorescence [5], where a beam of atoms 
interacts with a coherent laser field and all the electromagnetic modes of the vacuum [6]. 
Moreover, it is by now recognized that nearly all types of laser-atom interactions can be 
strongly affected by laser noise. Indeed one practical reason to this fact is the use, in 
experiments, of high powers which are obtained in pulsed operation, and thus at the expense 
of poorly stabilized laser beams. Furthermore, real atoms experience a fluctuating 
environment of many perturbing interactions and ideal lasers exist only in theoretical models, 
while the used laser sources are subjected to many types of fluctuations notably in phase, 
amplitude and frequency [7-10]. Other kinds of fluctuations due to collisional effects can 
affect the atomic transition frequencies [9,11]. Therefore, we cannot establish, without taking 
into account of the dissipative action of the environment and the statistical properties of the 
laser light, a rigorous comparison between theoretical predictions and experimental results. 

Different approaches to the dissipative dynamics of open quantum systems in strong 
external fields has been proposed and applied to the description of atoms under the influence 
of thermal noise [1-4]. While for the incoherence of laser field, a series of models, all based 
on 

so-called pre-Gaussian Markovian processes [7-10] has been used in order to describe the 
stochastic behavior of the laser field. It is important to mention a few technical features of this 
model. It is based on the two-state random telegraph. It is not a Gaussian model but rather 



"pre-Gaussian" and has a Gaussian limit [8]. Our choice of the random telegraph is based on 
the simplicity of this model that permits a unified treatment of different noisy laser in phase, 
amplitude and frequency. Several works have reported on the action of random process on a 
two-level system [7-10,12-15], particularly the evolution populations o nn and the ionization 
probability. 

In the present paper, we want to elucidate the role of quantum dissipation and laser 
noise on the atomic response. For this purpose, we derive a master equation, which provides a 
general framework for the dynamics of atoms interacting with strong laser noisy and thermal 
reservoir. 

The basic idea underlying the theoretical formalism of that paper is to take into 
account the exact dynamics of the interaction between atomic system and external field by 
employing the Floquet basis for the reduced system rather than the stationary unperturbed 
states 

[1,2]. The interaction of the system (laser-atom) with the reservoir will be treated by the time- 
dependent perturbation theory, this treatment leads to a generalized quantum master equation 
for the reduced density matrix. Such statistical fundamental equation firstly introduced in 
quantum optic by Burshtein [16-18], contains all information concerning, respectively, the 
atomic transition dynamics, the stochastic evolution of laser field fluctuations and the 
dissipative mechanisms. We are concerned here with an important theme of contemporary 
research, namely the interplay between quantum coherence and external noise. In fact, the 
destruction of quantum coherence by noise is central to many fields and is reflected in the 
large number of papers recently published on this subject [19-23]. 

Since our computations are made at exact resonance, where the effects of spontaneous 
emission are important [24] and for strong laser field, where the probability to realize an 



ionization of atoms is highest. Therefore we shall here extend our theory by the inclusion of 
the relaxation rates corresponding to the spontaneous emission and the ionization processes 
and present the corresponding numerical results. 

The theory is developed in section 2, by considering the case of strong laser-atom 
interaction in the presence of laser noise and dissipative effects, which are introduced by the 
coupling to the reservoir. Within the framework of the Floquet representation and the 
Markovian approximation, detailed theoretical calculations are feasible to obtain the reduced 
matrix density elements. The account of Floquet theory given here is rather brief, since the 
theory has been discussed at length in the recent literature (see e. g. refs. [25-29]). Moreover, 
the influence of ionization process on the response atomic function is presented. Numerical 
results concerning a model of two-level system is presented in section 3. At the end a 
summary of our results is given. 



II- THEORY 



We consider in this paper an atomic system, which interacts with an external classical 
laser field. Moreover, the system (laser-atom) is coupled to a quantified radiation field in 
thermal equilibrium. In the following we will consider the behavior of the atom coupled to a 
reservoir with many degrees of freedom. 

The aim of this paper is to provide a description of the dynamics in terms of the 
degrees of freedom of the atomic system alone by elimination of the reservoir variables. Since 
the atoms are driven strongly by an external laser field, our master equation is based on the 
atomic Floquet states rather than the unperturbed atomic states. 

The total Hamiltonian governing the dynamics of the coupled system of matter and 
radiation degrees of freedom takes the form 



where H A _ L (t) is the total Hamiltonian of the atomic system and the external laser field, 
without interaction with the reservoir, given by 



with H A denotes the stationary atomic Hamiltonian, V(t) the dipole interaction between the 
atomic system and the laser field and H SE the Hamiltonian of the simultaneous emission, 
which reads 



H(t) = H A _ L (t) + H I+ H R 



(1) 



H A _ L (t) = H A +V(t) + H 



SE 



(2) 



SE — 




(3) 



where T is a diagonal matrix composed by the Einstein coefficients of spontaneous 
emission process and defined by [29] 



^ LY S ^' (4) 

n'( n 



here y^, are the radiative decay rates. 

The Hamiltonian that describes the coupling between the matter degrees of freedom and the 
quantified radiation field may be written in the dipole approximation as 



H^jXy^+a+y;) (5) 



where y^ are the coupling constants. ajand aj are the quantum operators of creation and 

annihilation, z denotes the component of the dipole operator on the (OZ) axis. The free 
Hamiltonian of the reservoir is represented by 



H R =^£cOj 



f ♦ 1 A 
V L J 



(6) 



with (Oj is a frequency corresponding to the j eme mode of the free quantified radiation. 
The Von Neumann equation for the statistical operator p of the total system reads 



i*^ = [H(t),p] (7) 
at 



We introduce the interaction representation for treating the equation (7), and we put 



H (t) = H A _ L (t) + H R 



(8) 



that is considered as the time-dependent unperturbed Hamiltonian, the evolution operator 
corresponding to this Hamiltonian is given by 



U(t) = U (t)®U R (t) 



(9) 



and 



U(t) = 



eX p{-^> HA - L(t ' )dt '} 



exp(--H R t) 
n 



(10) 



where [. . .] is an operator of time ordering. In the interaction representation, the total density 
operator p(t) and the interaction Hamiltonian H x take the following form 



p(t) = U + (t)p(t)U(t) 



(ID 



and 



H^mo^ua) 



(12) 



and the dynamic equation (7 ) becomes 



i^=[H l( t),p(t)] 



at 



(13) 



We assume that the interaction between the atomic system and the reservoir is weak that the 
coupling constants y ■ — > and y?t = Cst for t — > oo . In these conditions, the equation (13) 

will be treated by the time-dependent perturbation theory. At the second order in Hp this last 
reads 

^ = Mu, (t),p(t )]~ f df [h : (t) , [h : (f ) , p(f )] ] (14) 

dt m h Jt <> 

In writing Eq. (14), it has been assumed that the interaction is adiabatically switched on at 
time (t — > - oo). Prior to this, the atomic system and the reservoir are uncorrelated and the 
total density matrix is given by the direct product 

p(t ) = o(t )®p R (15) 

where d(t ) is the reduced system density operator at initial time in the interaction 
representation and defined by the trace over the reservoir states. p R is the reservoir 
distribution function at equilibrium given by 



exp(-H R /K R T R ) .... 
P R = P * — (16) 



here Z R , T R are respectively the partition function and the temperature of the reservoir and K B 
is the Boltzman constant. 

We need to define the Hamiltonian of interaction between the atomic system and the laser 
field without its coupling with the bath. In the dipole approximation, it reads 



V(t) = e F zcos(tot + (p(t) + x(t)) 



(17) 



in the case of phase fluctuations, 



or 



V(t) = eF (l + x(t)) 



zcos 



(tot + (p(t)) 



(18) 



in the case of amplitude fluctuations, where CO is the laser frequency and e is the electron 
charge. F is the electric field amplitude (possibly fluctuating in magnitude) and <p(t) is the 

instantaneous phase of the laser (fluctuating around the mean value). In this paper, we use an 
intense laser field affected by a temporal stochastic process of jumps. These fluctuation 
mechanisms are described by the pre-Gaussian Markovian models [7-10]. In particular, we 
adopt the simplest example of two- state random telegraph, which is defined by x(t) = ± a, 
where a is the amount of the jump assigned to the stochastic signal. 

Since the telegraph process that we are considering here is Markovian, the 
conditioned probability density function associated with it, namely p(s, t |s , t ) , is shown 
to satisfy the following Chapman-Kolmogorov equation [7-10,30] 



dt 



p(s,t|s ,t ) = --p(s,t|s ,t ) + -p(-s,t|s ,t ) 



(19) 



here s is the initial state of random telegraph at the time t . 
In the compact form [31], the equation (19) write as 



(20) 



where W s s =^ 



-1 1 
1 -1 



is the relaxation matrix composed by the frequencies of telegraph 



jumps process, where s and s' are two different states of random telegraph (s = 1,2) and 
corresponding to the telegraph signal amplitude {-a, + a}. T denotes the dwell time (i.e., the 
mean time between interruptions) of the telegraph. In the following, and in presence of the 
noise, all physical operators, such reduced density operator, interaction Hamiltonian, dipole 
operator. . ., will be indiced by s in indication of the stochasticity influence. 
By elimination of the reservoir variables in Eq. (14), we have 



d.(t) = -^ fdt'Tr R [fi IS (t), [H IS (t'),o s (t')]] (21) 



After tracing on the reservoir variables, we combine the Chapman-Kolmogorov equation (20), 
which represents the stochastic evolution of the random telegraph, to the equation (21) 
representing the atomic dynamics. A master equation for the reduced density operator is 
derived in the interaction picture 

5 s (t) = £W s s d s ,(t) + D s (t) (22) 

s' 

where 

D s (t)=-^fdt'Tr R [H IS (t),[H IS (t'),o s (t')]] (23) 

h to 



is the time-dependent operator describing the dissipative effects induced by the coupling to 
the reservoir. In writing the expression of D s (t), we have replace up p s (t) ~5 s (t) ®p R to 
second order in the coupling constants by its zeroth-order approximate. Since H ls (t) is a 

periodic function in time, we explicitly construct the operator U (t). To this end, it is 
necessary to treat the interaction with the strong incoherent laser field exactly, and solve the 
corresponding Schrodinguer equation generated by the Hamiltonian H A _ L (t) by using the 
Floquet theory. According to Floquet's theorem [25-29], there exists a complete set of 
solutions labeled by quantum numbers a of the form 



i|C(t) =exp(-ie a t/ft) <L, s (t) (24) 



Where he a and <|> ajS (t)y are respectively the quasi-energies and the eigenstates of Floquet. 

The time-evolution operator U (t) for the matter degrees of freedom in the Floquet 
representation [1,2], is given by 

U (t) = exp( ie a t//z)|<L, s (t))(<L, s (0) | (25) 

In the interaction picture, the operators annihilation, creation and dipole respectively aj, 
a | and z take the form 



=U + (t)a J U(t) = exp(-ia) J t) &J , 
a/ = exp(i(Q J t)a+ 



(26) 
(27) 



and 



z s (t) = U + (t)zU(t) = £ exp(i8 ap (k)t)Z aPjS (k) | <L, (t))((() a , s (0) | (28) 



apk 



where 



a) — 

a M (k) = — | ) a) dtexp[-ikcot]((t) aiS (t)|z|(t) ps (t) 



(29) 



are the dipole matrix elements between Floquet states. The symbol 'prime' in the sum Eq. 
(28) indicates that only the triplets (oc,p\k) that verify the condition £ a -£ p +k/z(Q>0 will 
be considered, in the purpose to eliminate the degenerate Floquet eigenstates. In order to 
calculate the dissipation operator D s (t), we follow the methodology formulated in the 
Ref. [1], indeed, we have 



H IS (t) = ££ Yj 

a|3k j 



e ^ap.s a j +ti.C. 



(30) 



with H.c. is the hermitic conjugate, 



F„ M = *...(0) 4> M (0) 



(31) 



and 



£ap( k ) = ea -e p +kfico 



(32) 



where - °° < k < +oo . By taking into account of the (Eq. 30), D s (t) reads in the form 



6 s(t) = 11 Etff dt'^alaj^atOF^d.d')^, „ - A + (t,t')F pa , s F a , p , s a s (t')) 

apka'pv j 

+ (a J a;)(A + (t,t')F pa , s d s (t')F a , p , s - A(t,f)F ap , s F p , as d s (f))+H.c.} (33) 



with 



A(t,t') = Z aps (k)Z;, a , s (k')exp(i(E ap (k)-a) J )t)exp(i(8 a , p ,(k')-co j )t') (34) 



and 



) = N«D j ) = 



exp 



K B T R 



-1 



(35) 



N(C0j) is the photon number operator. Since the freedom degrees of the bath are infinite, we 

can make the substitution ^Yj--- — > JdcoJ(G))..., where J(to) is a function which is 

j 

proportional to the bath spectral density. In order to perform the integration, which is present 
in the expression of D s (t) Eq. (33), further conditions must be imposed on the reservoir in 
the purpose to prevent the energy, initially in the atomic system, from returning back from the 
heat bath to the system in any finite time, i.e., (treat the coupling of the reduced system to the 
reservoir as an irreversible process). At this stage, we make two approximations, 
i) Equation (33) contains o s (t') in the integral, and hence the behavior of the atomic system 

depends on its past history from t' = t to t'= t . The motion of the atomic system is 
however, damped by the coupling to the reservoir and damping destroys the knowledge of the 

past behavior of the system. Therefore the first assumption is that 5(t) depends only on it's 
present value c s (t) (Markovian approximation) [1,2]. 



ii) Let us consider an operator B of the bath and it's time correlation function ^B(t-t')B + ^ . 
Since the reservoir is assumed to be large and Markovian. Thus it is expected that 
^B(t-t')B + ^ will be nonzero for some time interval t-t' < x R , where x R is the correlation 

time of the reservoir. Interactions at times t and t' become progressively less correlated for 
t-t' » x R . The correlation function (^B(t-t')B + ^ is only maximum at t = t' . We can 

therefore tend the superior born of integration in Eq. (33) to the infinite ( t — > °° ). 

With these two approximations and by using the following expression for the initial 
time (t — > -oo), 

£df exp( i (co - e a . p . (k' )) f) = 2ti 5(a) - e a , p , (k' )) (36) 
The integro-differential equation reads 

D s (t) = 27l£ £j(8 a , p ,(k')){N(8 a , p ,(k'))(A'(t)F ap , s d s (t')F p a , s -A' + (t)F pa , s F a , p , s d s (t)) 

a|3k a'P'k' 

+ (l+N(e a , p ,(k')))(A' + (t)F pa , s o s (t)F a , p , s -A'(t)F ap>i F p , a , s 5 s (t))+H.c.} (37) 
Where 

A(t) = exp(i (e ap (k) -8 a , p ,(k))t )z ap , s (k) Z^k 1 ). (38) 
This last quantity is maximal for 



e ap( k )- e ap'( k/ ) = 2n ^ 



(39) 



Where n is a positive or negative integer, for the case of n = only terms such as 
(a, p\ k) = (a' , (3' , k' ) must be kept in Eq. (37). The equation (33) takes the final form 



D s (t) = ^n ap (k){N(8 ap (k))([F aM 5 s (t),F pas ] + [F aM ,5 s (t)F pa , s ]) 



a|3k 
+ 



(40) 



(l + N(8 ap (k)))([F pas d s (t),F aM ] + [F pa , s ,5 s (t)F aM ])} 



with 



Q ap (k) = 27Uj(8 ap (k))|Z aM (k) 



(41) 



By projecting on the Floquet basis j|(|) a >s (0))} , the master equation for the diagonal and the 
off- diagonal elements respectively S aa s (t) and o ap>s (t) read 



5««, s (t) = £w;d, (1/ , (t) + £ (M 7 „o m (t) - M ay a aa>s (t)) 



(42) 



and 



s«p,(t) = £w;o ap , s ,(t)-^ £(M aY +Mj 



5 ap>s (t) 



(43) 



where the coefficients M ap are defined by 



m «p =2 L{( 1 + N ( £ «p «)) a «P ( k ) + N ( e «P 00) a P« (k) } 



(44) 



and their solutions are given by 



S....(t) = E[«p(-Alt)L Ml S pp> , (0) 



(45) 



and 



where 



g aP, s (t) = L[exp(-A2t)] ay , ss ,d Yp ,, (0) 



y,s 



(46) 



Al a p jSS ' - W s 8 a p 



M pa -5 ap £M aT 



i J 



(47) 



and 



A2 ay, ss - = - W S S 8 ay + | 8 ay 8 SS , 



E (M an +M N ) 



(48) 



The theoretical expressions for populations and coherence of quasienergie states respectively 
Eq. (45) and Eq. (46) have to be transformed back into the atomic basis which yields G nn and 
a nn >. In Schrodinguer picture we then obtain 



o s (t) = - ^[h (t),o. (t)]+ £ w; o, (t) + D, (t) 



(49) 



with 



D s (t) = U D s U; 



(50) 



D s (t) = £a ap (k){(l+N(8 ap (k)))([R + ap , s (t),o s (t)R ap , s (t)]+ k p , s (t)a s (t),R ap , s (t)]) 

a|3k 

+ N(e ap (k)) ( [R ap , s (t), o,(t) R; p , s (t)]+ [R ap , s (t) o,(t), R + ap , s (t) ] ) } (51) 



where 



R 



a p fi (t) = U (t)F aM Uj(t) = exp(-i(e a -ep)t) |^ a , s (t))(^ p , s (t)| (52) 



The main difficulty of typical problems lies in the correct averaging of the matrix density 



over all realizations of noise. In fact, what is physically wanted is (o nn ), that is, the 



solution to the master equation in the atomic states and averaged over the ensemble of 



jumps of the implicit telegraph x (t). To obtain (o nn . ) one proceeds indirectly, by defining 



a marginal average o nn , s (t) , given by the equation 



(53) 



where g(s) is the initial probability distribution of the random process and o nn , s (t) the 
average value of a nn , (t) under the condition that x (t) is fixed at the value s at time t. 
By projecting on the atomic basis {| n)} , the master equation in the Schrodinguer picture then 
finally reads 



<W t )=E w ''<w( t )+E 



Dl m , n , s (t) + -H (t) m , n , s (t) 
n 



<w.(t)+ 



( i ^ 

LDl mm - « (t) — H n (t) mm , „ J mm . , (t) + V D2 „ (t)o m ,„, s (t) 
mm ,s v / j. U v ' mm ,s mm ,s v ' mmnn,sv / m n ,s v / 

m' y J m'n' 



(54) 



with the two terms responsible of the dissipation are defined by 



Di mn , s (t)=-^n aP (k){(i+N(8 aP (k))(t) ma , s (t) (t>: a , s (t)+ 

«Pk (55) 

N(e ap (k))(|> mM (t) (|>; M } 

and 

E>2 mm . 1I . 11 ..(t) = 2^Q ap (k){(l + N(e a p(k)))4 >m p..(t)<|>^ a ..(t)<|» 11 . a ..(t)<|»: |l ..(t) + 

apk (56) 

N(8 ap (k))(^ ma , s (t)(^ + m , p , s (t)(^ n , p , s (t)(^; a , s (t)} 

where (|) na s (t) are the Floquet states, which are projected on the atomic basis {| n}} . 

It is important to note that the general master equation (54) contains dissipative terms 
Eq. (55) and Eq. (56) that explicitly depend on time. This is the main difference to the usual 
optical Bloch equations. The physical interpretation to this fact is the strong distortion of the 
atomic dipole moment, which is induced by the external laser field. Since the atom couples to 
the environment via its dipole moment, the laser field also strongly influences the dissipation 
process [2]. 



Ill- Results and discussion. 



In this section we gather typical numerical results for the excitation and ionization of 
two-level atoms by strong laser fields in the presence of noise and dissipation mechanisms. 
To illustrate the effects of dissipation and laser noise on the atomic response, we present the 
evolution of atomic populations o nn (t) , which are obtained by numerical integration of the 
master equation (54). Our theoretical formalism is valid for the general case of multilevel 
systems but in order to keep the discussion simple; we will restrict our application to the 
two-level atoms for which a detailed study of the dissipative non-linear dynamics will be 
presented. A particular attention will be paid to the case of strong laser field, where the dipole 

operator is taken between the Floquet eigenstates { j^a.sCO) } rather than between 
unperturbed atomic states {I n}} . 




Having established the effects of strong laser noise on the atomic response and 
explored some features of different sources of noise (phase, amplitude and frequency). We 
concentrate our attention in this paper to the examination of quantum dissipation induced by 
coupling to the reservoir and when the noise is added to the laser field. We choose the inverse 
Rabi frequency Q. as time unit in the aim to analyze the obtained results in term of the noise 
magnitude. We are interested by a large light intensity such that the Rabi frequency is set to 
be equal to the atomic unit (Q. = 1 a.u.), This certainly is a very strong intensity. 




III. 1 Dissipative nonlinear quantum dynamics in the excitation of two-level systems. 



We begin by representing only the effects of quantum dissipation on the atomic 
response. Figures 1(a) and 1(b) show time evolution of populations of a two-level atomic 
system driven by strong coherent laser field sufficiently intense to remove a significant 
fraction of the population from the atomic ground state. One might think that the only 
consequence of a field this intense would be to lower the overall of the atom. The optical field 
is nearly resonant with allowed transition between discrete states of two-level atom. In 
absence of spontaneous emission decay and ionization effects, Figure 1(a) represents the 
atomic response without dissipation, the atomic system oscillates between the ground state 1 1) 
and some other discrete level | 2) and we have the ordinary Rabi oscillations. In Figure 1(b) 
the dissipation effect introduces a damping of Rabi oscillations. If damping effects are 
present, we expect that the Rabi oscillations will eventually become damped out and that the 
population's difference will approach some steady-state value for large scale of time. Hence, 
Rabi oscillations are not present in the steady state. 

It's interesting to note the presence of an irregular behavior on the oscillations of the 
two populations for a strong laser field; in fact we observe small oscillations which come to 
superpose to the Rabi oscillations, their amplitude is weak and disappears when the electric 
field strength F becomes small with respect to the atomic unit of field strength. These little 

oscillations represent the fast no-rotating variable phases exp(±i (u)+co nn , )t) in laser-atom 

interaction, which is treated in a non-perturbative way (Floquet theory). 

Figure 2, shows some typical results, in fact we take into account of both laser phase 
noise and dissipation influences. Taking a phase jump a = 0.4tc and three switching rates 
(Q.T = 0.1, 1 and 10). We display the time evolution populations on two column, in the 



column (A), only the effect of laser noise is considered, we remark in this case that for a 
switching rate (Q.T = 1), i.e., the noise frequency (1/T) is of the same size order that the Rabi 
frequency Q., a destruction of the atomic coherence is observed. The damping rate is strong 
and the relaxation to the steady state is rapid. The Rabi oscillations are restored when we 
consider the case of slow fluctuations (Q.T = 10) and fast fluctuations (CIT =0.1) and the 
damping rate is weak. The column (B) represents the same situation but by introducing 
moreover of the laser noise the dissipative terms (see Eqs (55) and (56)). A similar behavior is 
remarked as the column (A), but with damping rate more intense. For a switching rate 
(Q.T = 0.1), a partial destruction of the atomic coherence is induced by the dissipation effects. 
The kinetic of populations relaxation is more rapid that in the column (A). We remark that 
one of the effects of the quantum dissipation is the breaking of the atomic coherence 
especially for the case of (Q.T = 0.1) and the establishment of the dynamical localization 
regime for the case of (Q.T = 1). In Figure 3, we plot the time evolution of two-level atom 
populations in two situations, in first time, by neglecting the effect of quantum dissipation 
(column (A)) and considering only an amplitude telegraph noise and in second time we 
combine the two influences of noise and quantum dissipation. Taking an amplitude jump 
a = 0.1 a.u. and three different switching rates (Q.T = 0.1, 1 and 100). The column (A) shows 
pronounced quasi-oscillations, we remark a very weak damping at (Q.T = 1 and QT=100) and 
rapid relaxation for (Q.T = 0.1) with respect to the both cases of (Q.T = 1 and 100). In order to 
lead the system to the steady state we must use a large number of Rabi periods rather than in 
the case of phase fluctuations. The column (B) shows closely similar behavior that in column 
(A). The complicated structure of these curves is a consequence of the action of amplitude 
laser noisy on the reduced system dynamics. In fact, one observes a separation between the 
two occupation probabilities, every population G n (t) and o 22 (t) performs independently 



irregular oscillations, which converge to stationary state. In constraints to the case of phase, 
the addition of dissipation in Figure 3 (column (B)) introduces a weakness of the damping. 

Figure 4 illustrates the case of frequency fluctuations, this kind of noise is introduced 
by collisional effects, indeed the transition frequency CO21 can also fluctuate around its fixed 
value. The simplest model of such interruption collisions [9,11] assumes that the atomic 
transition frequency ©21 should be replaced by (O21 (t) = ©21 + x(t). By taking a jump 
parameter 

a = 0.1 a.u. and three different frequency switching rates (Q.T = 1, 10 and 100). We remark in 
the column (A) damped quasiperiodic oscillations. The case of (Q.T = 10) corresponds to 
strong damping without any convergence to a steady state. While the relaxation to an 
equilibrium state of value 1/2 is clear for a switching rate (Q.T = 1). The damping becomes 
weak for (Q.T = 100) and two independent beats phenomena are observed. The complicated 
time evolution of populations is a result of Rabi oscillation interference. In the column (B) 
where we take into account of quantum dissipation, the two populations relax to equilibrium 
state. The thermal noise induced by coupling to the bath introduces a complete destruction of 
the atomic coherence. The kinetic of relaxation and rate damping decrease from the case of 
fast fluctuations (Q.T = 0.1) to the slow fluctuations (Q.T = 100). 

In Figures 2 and 4, the comparison between the two columns (A) and (B) shows that 
the dissipation which behaves as a noise (thermal noise) leads the system to an equilibrium 
state with rapid kinetic of relaxation. The damping rates become large when we introduce the 
dissipation terms (column B). In other hand, we clearly see the destruction of atomic 
coherence, which increases, when we take into account of dissipation . The dynamical 
localization regime appears for phase and frequency noises. However Figure 3 shows 
important asymmetries. This behavior is justified by the fact that in the case of amplitude 



fluctuations, the jump parameter a, assigned to stochastic process, appear in term of laser 
intensity F (1 ± a ), while in the case of phase noisy , the dependence occurs in term of 

exp(±ia). 



III. 2 Dissipative nonlinear quantum dynamics in the ionization of two-level systems. 

Having established a formal framework for the excitation of atoms by laser noisy in 
presence of the reservoir action and explored some of its general predictions. We turn now to 
the examination of ionization effects on the populations and the illustration of the 
modifications generated by the different kinds of noise and quantum dissipation on the 
ionization rates. In order to analyze the ionization effects we adopt the extended two-level 
system model proposed by Yeh and Eberly [14,32]. The computation of ionization probability 
is made by the incorporation of responsible term of ionization [-x EC o m s , with 

x ec = R mc &mn + 1/ 2(R mc + R nc )(l - ) and R nc is the relaxation rate from the excited state 
|n> to the continuum |c>] in the motion equation (54). The trace of a over a complete set 



atomic states leads to the expression P ion (t) = 1 - >, <5 nn for the total ionization probability of 



the system [4,13,33]. 

We begin by showing successively the effects of laser noise and reservoir dissipation 
on the ionization probability. As illustrated in Figure 5, we have plotted total ionization as a 
function of Rabi periods. Taking an intense electric laser field such as (Q. = 1 a.u.) and a 
resonant laser frequency. The four curves of Figure 5 correspond respectively to the situations 
where: noise and dissipation are neglected, only the noise is considered, only dissipation 
effect is retained and both noise and dissipation exist. In Figure 5(a) we have considered a 



2 




phase noise of a = OAn, the results depend on the fluctuations time scale (1/T) compared to 
the other characteristic time scales of the problem such the Rabi frequency Q.. The minimum 
variations of the ionization probability are obtained, where we neglect both noise and 
dissipation. When we take into account of dissipation, the ionization probability increases. 
For large number of Rabi period, we remark that the phase fluctuations effect is more 
important that the reservoir action, in fact, the ionization probability rapidly increases when 
we introduce phase fluctuations corresponding to the case of (Q.T = 0.1). In presence of 
dissipation terms, the ionization probability variations remain closely constant. We conclude 
that noise and dissipation rapidly leads the atom to the ionization states. Figure 5(b) shows the 
case of amplitude noise, where the jump parameter is taken a = 0. 1 a.u. and a switching rate 
(Q.T = 1). A close behavior is observed as Figure 5(a) but the amplitude noise effect is very 
weak. The technique that we have used will be applied subsequently to frequency noise. 

A plot that gives a pictorial sense of how ionization proceeds in time is given in 
Figure 6. The total probability of ionization P ion (t) as a function of the Rabi period is plotted 
together with the occupation probabilities of bound states. The oscillations in these curves 
reflect the Rabi oscillations of the atom between the resonantly coupled states |l) and |2). 

These oscillations are damped by ionization in a few Qt's. This behavior is well known from 
the study of bound states coupled by an intense field. The column (A) of Figure 6 shows the 
response of two-level in presence of phase noise and (B column) by considering moreover the 
effect of dissipation. The same parameters are taken as Figure 2. The populations, which have 
not been lost through direct ionization to the atomic continuum, oscillate in the same manner 
that in absence of ionization effect, but there is a progressive decay to a zero probability. The 
ionization probability can be viewed as dominant in a few Rabi periods and increases rapidly 



in time when the quantum dissipative effects are considered. Figure 7 displays the same 

behavior that Figure 6 but for an amplitude noise. 

CONCLUSION 

In this paper we have investigated at length the non-linear dynamics of dissipative 
quantum atomic systems subjected to the action of heat bath and periodically driven by strong 
laser field which is affected by classical noise. We have derived and solved a master equation 
for atoms in strong noisy laser fields and in presence of reservoir dissipative effects. Such 
equation based on the Floquet states rather than the unperturbed atomic states has given 
typical and interesting results concerning the atomic dynamics. In fact, we have demonstrated 
how the master equation formalism, the Floquet theory, the pre-Gaussian models of laser 
noise and the Markovian coupling of quantum system to an environment can be combined 
together in order to tackle a general theoretical formalism and a powerful tool for the detailed 
analysis of the interaction of an atomic system with intense incoherent laser field and with a 
large reservoir. 

We have shown that the dissipation terms, which are time dependent respect to those 
in the Bloch equations, force the system to settle to some " preferred states " it is the 
dynamical localization regime observed in the cases of phase and frequency noises, as we 
have explored in this paper. Moreover under the action of these decay mechanisms, the 
atomic system exhibits different regimes such as the destruction of coherence, the relaxation 
to equilibrium state. In general, the strength of damping and the kinetic of relaxation increase 
with the addition of dissipation effects but the case of amplitude fluctuations show a different 
behavior. We have also analyzed the modifications induced by ionization effects. 
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Figure captions 



Figure 1. Populations o nn versus time (in units of inverse Rabi frequency Q.) for 
two-level atoms, resonantly excited by an intense electric laser field, such that the Rabi 
frequency is set Q. = 1 a.u. The emission spontaneous coefficient is y 2 \ = Q./10 5 a.u. 

(a) Both effects of noise and dissipation are neglected (coherent laser and no coupling 
to the reservoir). 

(b) In absence of noise and the dissipation is considered. 

Figure 2. Populations o nn versus time (in units of inverse Rabi frequency Q.) for 
two-level atom resonantly excited by random telegraph phase noise, successive frames are for 
different values of phase switching rates Q.T = 0.1, 1 and 10. We use strong laser field, such 
that the Rabi frequency is set Q. = 1 a.u. The phase jump parameter is a = 0.471, the emission 
spontaneous coefficient is y 2 i = £2/10 5 au. 

- The column (A) represents the effects of phase noise. 

- The column (B) represents the same situation as column (A) but by adding the dissipative 
effects. 

Figure 3. Same as Figure 2., but for an amplitude noise with the jump parameter is 
a = 0.1 a.u. and three switching rates Q.T = 0.1, 1 and 100. 

Figure 4. Same as Figure 3., but for a frequency noise. 



Figure 5. Ionisation probability P ion (t) versus time (in units of inverse Rabi frequency 
Q.). Same parameters as the previous figures are used. The The relaxation rate from bound 
states to the continuum is R 2c = Q/100 . 

(a) Dotted line: P ion (t) in absense of noise and dissipation. Dashed line: only the effect of 

phase noise for Q.T = 0.1. Dashed-Dotted line: only the effect of dissipation is considered and 
Solid line: both effects of noise and dissipation are considered. 

(b) Same as (a), but for an amplitude noise with Q.T = 1. 

Figure 6. The columns A and B are, respectively, same as Figure 2(A) and 
Figure 2(B), but take into account of the ionisation process represented by the probability 
P ion (t). The relaxation rate from bound states to the continuum is R 2c = Q/100 . 

Figure 7. The columns A and B are, respectively, same as Figure 3(A) and 
Figure 3(B), but take into account of the ionisation process represented by the probability 
Pion(t). The relaxation rate from bound states to the continuum is R 2c = Q/100 . 
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